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I.  INTRODUCTION 


The  problem  of  electromagnetic  scattering  from  a  perfectly  conducting 
surface  S  has  been  solved  by  first  placing  an  unknown  electric  current  J 
on  this  surface,  next  writing  the  electric  field  integral  equation  for  J, 
and  finally  numerically  solving  this  equation  by  the  method  of  moments  [1], 
[2].  The  solution  thus  obtained  will  be  accurate  if  S  encloses  neither  a 
resonant  cavity  nor  a  cavity  with  an  extremely  small  aperture,  and  if  the 
maximum  dimension  of  S  is  not  more  than  a  few  wavelengths.  The  surface  S 
may  be  two-dimensional  or  three-dimensional.  Of  course,  if  S  is  two- 
dimensional,  then  it  is  infinitely  long  in  one  direction,  and  we  require 
only  that  its  contour  in  the  plane  transverse  to  this  direction  be  not 
more  than  a  few  wavelengths  long. 

If  S  forms  a  cavity  with  a  very  small  aperture,  then  the  above 
method  will  fail  to  accurately  determine  the  field  inside  this  cavity 
because  this  field,  being  the  sum  of  the  incident  field  and  the  field  due 
to  J,  will  be  very  small.  The  field  due  to  J_  nearly  cancels  the  incident 
field.  Thus,  a  small  percentage  error  in  £  will  give  rise  to  a  large  per¬ 
centage  error  in  the  field  inside  the  cavity. 

The  field  inside  a  cavity  with  a  small  aperture  has  been  accurately 
determined  by  expressing  it  as  the  sum  of  the  short-circuit  field  and  the 
field  due  to  an  electric  current  source  which  is  the  negative  of  the  short- 
circuit  current  on  the  shorted  aperture  [3] .  This  electric  current  source 
radiates  in  the  presence  of  the  conducting  surface  S.  The  short-circuit 
field  is  the  field  that  would  exist  if  the  aperture  were  closed  by  a 
perfect  conductor.  The  short-circuit  current  is  the  electric  current  that 
would  exist  on  S  if  the  aperture  were  closed  by  a  perfect  conductor. 


An  alternative  method  called  the  generalized  network  formulation 
for  aperture  problems  [4] -[6]  has  been  proposed  for  the  case  in  which  S 
forms  a  cavity  with  a  small  aperture.  In  this  method,  the  aperture  is 
closed  with  a  perfect  conductor,  a  sheet  of  unknown  magnetic  current  M 
is  placed  on  one  side  of  the  shorted  aperture,  and  -M  is  placed  on  the 
other  side.  The  resulting  situation  will  be  the  same  as  that  of  the 
original  scattering  problem  if  the  component  of  magnetic  field  tangent 
to  the  shorted  aperture  is  continuous  across  the  shorted  aperture. 

Following  the  method  of  moments,  a  linear  combination  of  expansion  func¬ 
tions  is  substituted  for  M  into  the  equation  that  expresses  continuity  of 
the  tangential  magnetic  field  across  the  shorted  aperture  and  then  this 
equation  is  solved  numerically  for  the  coefficients  of  the  expansion  func¬ 
tions.  The  generalized  network  formulation  for  aperture  problems  is  easily 
implemented  for  an  aperture  in  an  infinite  conducting  plane  because  the 
electric  current  induced  on  the  complete  conducting  plane  by  each  expansion 
function  for  M  can  be  accounted  for  by  imaging  that  expansion  function  about 
the  plane  [5],  [6].  The  complete  conducting  plane  is  the  plane  with  the 
aperture  shorted. 

When  the  aperture  is  in  a  curved  surface,  the  generalized  network 
formulation  for  aperture  problems  is  not  as  easy  to  implement  because  it  is 
difficult  to  obtain  the  electromagnetic  field  produced  by  each  expansion 
function  for  M  placed  on  either  side  of  the  shorted  aperture  and  radiating 
in  the  presence  of  the  complete  conducting  surface.  The  complete  conduct¬ 
ing  surface  is  the  conducting  surface  S  with  its  aperture  shorted.  The 
problem  of  obtaining  the  electromagnetic  field  due  to  such  an  expansion 
function  so  radiating  is  a  scattering  problem  in  which  the  impressed 
source  is  on  the  scatterer.  Two  static  problems  of  this  type  were  con- 


sidered  in  [7].  In  the  first  problem,  the  objective  was  to  find  the 
electrostatic  field  due  to  an  electrostatic  dipole  on  the  conducting 
surface  of  a  sphere.  The  dipole  was  normal  to  the  sphere.  In  the 
second  problem,  the  objective  was  to  find  the  magnetostatic  field  due 
to  a  magnetostatic  dipole  tangentially  placed  on  the  conducting  sur¬ 
face  of  a  sphere.  In  [7],  each  of  these  two  static  scattering  problems 
was  solved  by  a  method  in  which  a  pseudo-image  was  used.  We  call  this 
method  the  pseudo-image  method. 

The  problem  of  electromagnetic  scattering  from  an  infinitely  long, 
infinitesimally  thin,  perfectly  conducting  cylindrical  surface  with  an 
infinitely  long  but  narrow  slot  is  stated  in  Section  II.  Solution  by 
means  of  the  generalized  network  formulation  for  aperture  problems  is 
outlined  in  Section  III.  In  Sections  II  and  III,  the  electromagnetic 
excitation  is  from  impressed  sources.  The  field  due  to  the  impressed 
sources  radiating  in  the  presence  of  the  complete  conducting  surface  is 
called  the  short-circuit  field  due  to  the  impressed  sources.  This  field 
is  evaluated  in  Section  IV.  Only  one  expansion  function  for  M  was  used 
in  Section  III.  In  Section  V,  the  pseudo-image  method  is  used  to  obtain 
the  short-circuit  field  due  to  this  expansion  function  placed  on  either 
side  of  the  shorted  aperture.  In  Section  VI,  coefficients  appearing  in 
previously  derived  expressions  for  the  short-circuit  fields  due  to  the 
impressed  sources  and  to  the  expansion  function  for  M  are  evaluated  for 
TM  plane  wave  excitation.  In  Section  VII,  the  coefficient  of  the  expan¬ 
sion  function  for  M  is  evaluated  for  TM  plane  wave  excitation. 


STATEMENT  OF  THE  PROBLEM 


The  combination  of  a  known  impressed  electric  current  source 

and  a  known  impressed  magnetic  current  source  Mlm^  radiates  in  a  homogeneous 

medium  with  permeability  U  and  permittivity  £  to  produce  an  incident  electro- 
inc  inc 

magnetic  field  (E  ,  H  )  that  is  independent  of  z.  An  infinitely  long 
perfectly  conducting  cylindrical  surface  called  S  is  now  placed  in  the 
medium.  The  axis  of  S  is  the  z  axis  so  that  the  electromagnetic  field 
radiated  by  Mimp)  in  the  presence  of  S  is  independent  of  z.  This 

field  is  called  (E,  H) . 

Our  purpose  is  to  determine  (E,  H)  when  S  is  infinitely  thin  and,  as 
viewed  in  the  xy  plane,  appears  as  the  contour  shown  in  Fig.  1.  This  con¬ 
tour  is  a  chain  of  N-l  straight  line  segments.  The  jth  segment  extends  from 
the  point  to  the  point  The  (x,y)  coordinates  of  the  points  P^  and 

PN  are  (0,w)  and  (0,-w),  respectively.  The  plane  strip  of  area  whose  pro¬ 
jection  on  the  xy  plane  is  the  straight  line  segment  between  points  P^ 
and  P^  in  Fig.  1  is  called  the  aperture.  If  there  were  no  gap  between  the 
points  P^  and  P^,  then  S  would  enclose  a  cavity.  If  the  gap  were  small,  the 
straightforward  method  via  solution  of  the  electric  field  integral  equation 
for  the  induced  electric  current  on  S  would  fail  to  accurately  determine 
(E,  H)  inside  the  cavity .  This  method  would  fail  because  it  gives  (E,  H)  as 
the  sum  of  two  fields  that  nearly  cancel  each  other  inside  the  cavity, 

(Einc,  HinC)  and  the  field  due  to  the  induced  electric  current  on  S.  The 
method  of  solution  described  in  the  next  section  accurately  determines 


(E,  H)  when  the  gap  is  small. 


APERTURE 


The  original  problem  as  viewed  in  the  xv  plane.  The  two- 
dimensional  sources  I10*3  and  Mimp  radiate  in  the  presence 
of  the  cylindrical  surface  S. 


The  original  problem  of  Fig.  1  is  replaced  by  the  equivalent  problem 
of  Fig.  2.  The  situation  of  Fig.  2  was  obtained  by  closing  the  aperture  with 
a  perfectly  conducting  plane  strip,  placing  the  magnetic  current  sheet  M  on 
the  left-hand  side  of  the  strip,  and  placing  -M  on  the  right-hand  side  of 
the  strip.  Here, 

M  =  u  x  E  (1) 

—  — x  — 

where  u  is  the  unit  vector  in  the  x  direction  and  E  is  the  electric  field 
— x  — 

in  the  aperture  in  Fig.  1.  The  combination  of  S  and  the  plane  strip  that 

SC  sc 

short  circuits  the  aperture  is  called  S  in  Fig.  2.  The  surface  S  divides 
space  into  two  regions,  an  external  region  called  region  a  and  an  interior 
or  cavity  region  called  region  b. 

If  there  was  no  electric  current  on  the  shorting  strip  in  Fig.  2, 
then  this  strip  could  be  removed.  After  removal  of  the  shorting  strip, 
the  magnetic  current  sheets  M  and  -M  would  cancel  each  other  resulting  in 
the  situation  of  Fig.  1.  Therefore,  the  situation  of  Fig.  2  will  be 
equivalent  to  that  of  Fig.  1  if  there  is  no  electric  current  on  the  short¬ 
ing  strip  in  Fig.  2.  There  will  be  no  electric  current  on  this  shorting 
strip  if  the  tangential  magnetic  field  immediately  to  the  left  of  the 
shorting  strip  is  equal  to  the  tangential  magnetic  field  immediately  to 
the  right  of  the  shorting  strip.  Hence,  we  write 

Ha  (Jimp,  Mimp)  +  Ha  (0,M)  =  H^  (0  ,-M)  on  A  (2) 

—tan  —  —  —tan  —  —tan  — 

where  A  is  the  surface  of  the  shorting  strip.  The  subscript  tan  in  (2) 
denotes  the  component  tangent  to  A.  In  (2),  M^5)  is  the  magnetic 


field  in  region  a  due  to  (Jimp,  Mimp)  radiating  in  the  presence  of  SSC, 

H  (0,  M)  is  the  magnetic  field  in  region  a  due  to  M  placed  on  the  left- 

hand  side  of  the  shorting  strip  and  radiating  in  the  presence  of  S  ,  and 
b  . 

H  (0,  -M)  is  the  magnetic  field  in  region  b  due  to  -M  placed  on  the  right- 
hand  side  of  the  shorting  strip  and  radiating  in  the  presence  of  S  .  Each 
magnetic  field  in  (2)  has  two  arguments.  The  first  argument  is  an  electric 
current,  and  the  second  argument  is  a  magnetic  current.  Actually,  the  left- 
hand  side  of  (2)  should  be  evaluated  not  exactly  on  A,  but  immediately  to 
the  left  of  A  in  region  a.  Moreover,  the  right-hand  side  of  (2)  should  be 
evaluated  not  exactly  on  A,  but  immediately  to  the  right  of  A  in  region  b. 
Equation  (2)  ensures  that  the  situation  of  Fig.  2  is  equivalent  to  that  of 


Fig.  1, 


.  * . 


Since  the  situation  of  Fig.  2  is  equivalent  to  that  of  Fig.  1,  the 


electromagnetic  field  in  Fig.  2  is  equal  to  the  unknown  electromagnetic 
field  (_E,  H)  in  Fig.  1.  The  field  (j£,  H)  in  Fig.  2  is  partitioned  into 
(Ea,  Ha)  and  (Eb,  Hb)  where  (E3,  H3)  is  the  field  in  region  a,  and  (Eb,  Hb) 


is  the  field  in  region  b. 


(E,  H)  = 


ci  a 

(E  ,  H  ),  region  a 


(Eb,  Hb),  region  b 


E3  =  E3(Jimp,  Mlmp)  +  Ea(0,  M) 
H3  =  Ha(Jlmp,  Mlmp)  +  Ha(0,  M) 


Eb  =  Eb(0,  -M) 


Hb  -  Hb(0,  -M) 


The  magnetic  fields  on  the  right-hand  sides  of  (5)  and  (7)  were  defined 
in  the  previous  paragraph.  The  E's  on  the  right-hand  sides  of  (4)  and 
(6)  are  the  corresponding  electric  fields.  Namely,  Ea ( ,  M^mp)  is 
the  electric  field  in  region  a  due  to  (Jimp,  Minip)  radiating  in  the  pre- 

SC  3. 

sence  of  S  ,  _E  (0,  M)  is  the  electric  field  in  region  a  due  to  M  placed 

on  the  left-hand  side  of  the  shorting  strip  and  radiating  in  the  presence 
sc  b 

of  S  ,  and  ji  (0,  -M)  is  the  electric  field  in  region  b  due  to  -M  placed 
on  the  right-hand  side  of  the  shorting  strip  and  radiating  in  the  presence 

SC 

of  S  .  Each  electric  field  on  the  right-hand  sides  of  (4)  and  (6)  has 
two  arguments.  The  first  argument  is  an  electric  current,  and  the  second 
argument  is  a  magnetic  current. 

Some  fields  on  the  right-hand  sides  of  (4) — (7 )  depend  on  M.  To 
determine  M,  we  appeal  to  (2) .  Rearranging  the  terms  in  (2) ,  we  obtain 

-Ha  (0,  M)  -  H*  (0,  M)  =  Ha  (Jimp,  Mimp)  on  A  (8) 

—tan  —  —tan  —  —tan  —  — 

Next,  M  is  approximated  by 

M  =  VMX  (9) 

where  V  is  an  unknown  complex  constant  to  be  evaluated  later  and  is 

an  expansion  function  that  will  be  specified  later  when  the  incident 

j.nc  inc  _  _ 

field  (E  ,  H  )  is  specialized  to  a  TM  field.  The  expansion  function 

will  be  tangent  to  SSC  and  independent  of  z  where  z  is  the  coordinate 

measured  perpendicular  to  the  xy  plane.  The  magnetic  fields  in  (8)  will 

also  be  independent  of  z  so  that  it  will  suffice  to  satisfy  (8)  only  on 

the  projection  of  A  in  the  xy  plane.  This  projection  is  the  straight  line 

segment  between  the  points  P^  and  in  Fig.  2.  Substituting  (9)  into  (8), 

taking  the  scalar  product  of  (8)  with  K.  ,  and  then  integrating  this  product 


from  -w  to  w  with  respect  to  y  at  x=0,  we  obtain 


(Ya  +  Yb)V  =  I 


(10) 


where 


w 


Ml  *  Ha(0,  M^)dy 


(ID 


-w 


w 


Mi  •  H  (0,  M^dy 


(12) 


-w 

w 


Ml  •  Ha(Jlmp,  Mimp)dy 


(13) 


-w 


Actually,  the  magnetic  fields  in  (11)  and  (13)  are  evaluated  not  exactly 

at  x=0  on  the  shorting  strip,  but  immediately  to  the  left  of  the  shorting 

strip  in  region  a.  Moreover,  the  magnetic  field  in  (12)  is  evaluated  not 

exactly  at  x=0  on  the  shorting  strip,  but  immediately  to  the  right  of  the 

shorting  strip  in  region  b. 

3.  b 

If  Y  ,  Y  ,  and  I  can  be  determined,  then  (10)  can  be  solved  for  V. 
Substitution  of  (9)  into  (4)-(7)  gives 


E 

H 

E 


a 


a 

b 


H 


=  Ea(Jlmp,  Mimp 

-  Ha(Jimp,  Mimp 

-  -  VEb(0,  M:) 

=*  -  VHb(0,  Mx) 


)  +  VE  (0,  Mx) 


)  +  VH  (0,  Mx) 


(14) 

(15) 

(16) 
(17) 


The  electromagnetic  fields  on  the  right-hand  sides  of  (14)-(17)  are 

obtained  in  the  next  two  sections.  The  magnetic  fields  so  obtained  are 

3  b 

then  substituted  into  expressions  (11) -(13)  for  Y  ,  Y  ,  and  I. 
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IV.  THE  SHORT-CIRCUIT  FIELD  DUE  TO  THE  IMPRESSED  SOURCES 


In  this  section,  the  electromagnetic  field  (JSa(J*m*>,  M*0^) , 
Ha(J1InP,  which  appears  in  (14)  and  (15)  is  evaluated.  This  field 

is  the  field  in  region  a  due  to  the  impressed  sources  (J^™^,  M10^)  radi- 

SC 

ating  in  the  presence  of  S  .  The  xy  plane  view  of  the  situation  in 
which  this  field  exists  is  shown  in  Fig.  3.  It  is  evident  from  Fig.  3 
that 


Ea(J1"!>,  M1”*’)  -  EinC 

+  E(JinC, 

0) 

in  region  a 

(18) 

Baalmp,  Mi,np)  -  H1”0 

+  H(JlnC, 

0) 

in  region  a 

(19) 

where  £lnc  is  the  electric  current  induced  on  Ssc  by  (Jlm*>,  Mim^) .  In 

(18)  and  (19),  (Einc,  H^nc)  is  the  electromagnetic  field  due  to  Mirap) 

radiating  in  the  presence  of  the  homogeneous  medium  characterized  by  (u»e), 

and  (E(JinC,  0),  H(JinC,  0))  is  the  electromagnetic  field  due  to  the  elec- 
Anc  _ 

trie  current  radiating  in  the  presence  of  the  same  homogeneous  medium. 

Although  the  right-hand  sides  of  (18)  and  (19)  are  defined  in  all  space, 
the  left-hand  sides  of  (18)  and  (19)  are  defined  only  in  region  a  so  that 
(18)  and  (19)  can  be  valid  only  in  region  a. 

Since  S  is  perfectly  conducting  in  Fig.  3, 


£a  (JimP>  MimP) 
—tan  —  — 


on  S 


sc 


(20) 


where  the  subscript  tan  denotes  the  component  tangent  to  S 


sc 


Because 


the  field  in  (20)  is  given  by  (18),  we  obtain  the  following  equation 

,  Tinc 
for  J 


-E  Uinc,  0)  -  e*™; 
—tan  —  —tan 


sc 


on  S 


(21) 


The  electric  current  in  (21)  is  approximated  by 


jilic 


N  ^ 

I  Iinc  J . 

j«l  3  -3 


(22) 


where  {l*nc}  are  unknown  constants  and  LJ }  are  expansion  functions  that 


will  be  specified  later  when  the  incident  field  (El^nc,  H^nc)  is  specialized 


to  a  TM  field.  Each  of  the  expansion  functions  { J  }  will  be  on  ,  tangent 


.sc 


to  S  ,  and  independent  of  z. 

Equation  (22)  is  substituted  into  (21)  and  then  the  scalar  product 


—  v  SC 

of  (21)  with  is  integrated  around  the  contour  of  S  in  the  xy  plane. 


Letting  i  »  1,2,...,  N,  we  obtain  N  simultaneous  equations  which  are  written 
in  matrix  form  as 

-W-i  nn 

(23) 


z^inc  _  +inc 


where  Z  is  a  square  matrix  of  order  N  and  each  of  the  quantities  I 


inc 


and 


-►inc  -►inc 

V  is  a  column  vector  of  N  elements.  The  jth  element  of  I  is  the 


.inc 


constant  1^  that  appears  in  (22).  The  ijth  element  of  Z  is  called  Z „  and 


is  given  by 


N+l 


'ij 


1±  *  I(ij»  0)dt 


(24) 


where  the  variable  of  integration  t  is  the  arc  length  along  the  contour 
of  S  in  the  xy  plane.  Starting  from  zero  at  the  point  in  Fig.  3,  the 

variable  t  increases  as  the  contour  is  transversed  in  the  clockwise  direc¬ 
tion,  attaining  the  value  tj^+1  upon  arrival  back  at  the  starting  point  P^. 
The  ith  element  of  V*nc  in  (23)  is  called  V*nc  and  is  given  by 


After  the  elements  of  Z  and  V  have  been  obtained,  (23)  can  be 


-*inc  inc 

solved  for  I  .  The  electric  current  J  will  then  be  given  by  (22) 


in  which  I  is  the  jth  element  of  1  .  Substitution  of  (22)  into 


(18)  and  (19)  gives 


Ea(Jimp,  Mimp)  -  Einc  +  l  lj-nc  E(J.,  0)  in  region  a  ( 

j-1  2  2 


N 


rinc 


and 


*  H^nC  +  £  I'!'110  HOJ.,  0)  in  region 

1-1  3  3 


H  (J 


Substituting  (27)  into  (13),  we  obtain 

N 


w 

r 


M  •  Hlncdy  +  l  I*nc  M..  •  H~(J.,  0)  dy 
“1  -  4”i  J  J-1  -  -J 


„inc 


-w 


j-1 


-w 


where  the  superscript  denotes  magnetic  field  evaluation  at  x-0  , 
that  is,  immediately  to  the  left  of  the  shorting  strip. 


V.  THE  PSEUDO-IMAGE  METHOD  OF  OBTAINING  THE  SHORT-CIRCUIT  FIELDS  DUE 


TO  THE  MAGNETIC  CURRENT  EXPANSION  FUNCTION 


In  this  section,  the  electromagnetic  fields  (Ea(0,  M^),  Ha(0,  M^)) 


and  (E^(0,  M^),  H^(0,  M^))  of  (14)-(17)  are  evaluated  by  means  of  a  method 


called  the  pseudo-image  method.  Here,  E  (0,  M^)  and  H  (0,  M^)  are,  respec 
tively,  the  electric  and  magnetic  fields  in  region  a  due  to  the  magnetic 


current  sheet  placed  on  the  left-hand  side  of  the  shorted  aperture  and 


radiating  in  the  nresence  of  S3C.  Moreover,  ^  (0,  M^)  and  H^(0,  M^)  are. 


respectively,  the  electric  and  magnetic  fields  in  region  b  due  to  the 


magnetic  current  sheet  placed  on  the  right-hand  side  of  the  shorted 

SC 

aperture  and  radiating  in  the  presence  of  S 

Figure  4  shows  the  xy  plane  view  of  the  situation  in  which  (Ea(0,  M^), 
Ha(0,  M^))  exists.  Since  the  perfectly  conducting  surface  SSC  completely 
isolates  region  b  from  region  a,  any  source  can  be  placed  in  region  b 
without  affecting  the  electromagnetic  field  (Ea(0,  M^),  Ha(0,  M^))  in 
region  a.  We  place  the  magnetic  current  sheet  on  the  right-hand  side 
of  the  shorted  aperture  as  shown  in  Fig.  5. 

The  field  in  region  a  of  Fig.  5  is  the  same  as  that  in  region  a 

3  3 

of  Fig.  4,  namely,  (E  (0,  M^) ,  ]i  (0,  M^)).  The  magnetic  current  sheet  on 
the  right-hand  side  of  the  shorted  aperture  in  Fig.  5  is  called  the  pseudo¬ 
image  of  the  magnetic  current  sheet  on  the  left-hand  side  of  the  shorted 
aperture.  The  pseudo-image  annihilates  the  primary  electric  field  exci¬ 
tation  on  the  shorted  aperture.  This  annihilation  should  considerably 
reduce  the  electric  current  induced  on  S  and,  hopefully,  improve  the 

behavior  of  this  current  so  that  it  can  be  more  accurately  approximated  by 
a  l:\iear  combination  of  simple  expansion  functions. 

It  is  evident  from  Fig.  5  that 

Ea(0,  -  E(0,  21^)  +  E(Ja,  0)  in  region  a  (29) 

and 

Ha(0,  Ml)  -  H(0,  2Mj_)  +  H(Ja,  0)  in  region  a  (30) 

3  SC 

where  J_  is  the  electric  current  induced  on  S  in  Fig.  5,  and  2M^  repre¬ 
sents  the  combination  of  the  two  magnetic  current  sheets  in  Fig.  5,  each 
of  density  M^.  On  the  right-hand  sides  of  (29)  and  (30),  (E(Ja,  0),  H(J  ,  0)) 
is  the  electromagnetic  field  due  to  the  electric  current  ^J3,  and  (E(0,  2M^)» 


y  REGION  a 


Fig.  5.  The  xy  plane  view  of  the  equivalent  situation  in  which  the 
electromagnetic  field  Ha(0,H^))  exists.  The 

magnetic  current  sheet  immediately  to  the  right  of  the 
shorted  aperture  is  called  the  pseudo-image  of  the  one 
immediately  to  the  left. 


H(0,  2M^))  is  the  electromagnetic  field  due  to  the  magnetic  current  2H^. 
In  the  previous  sentence,  both  J_  and  2M^  radiate  in  the  presence  of  the 
homogeneous  medium  characterized  by  (y,e).  The  right-hand  sides  of  (29) 
and  (30)  are  the  fields  everywhere  in  Fig.  5,  whereas  the  left-hand  sides 
of  (29)  and  (30)  are  defined  only  in  region  a  so  that  (29)  and  (30)  can 
be  valid  only  in  region  a. 

3  SL 

In  order  to  determine  _E  (0,  M^)  and  H  (0,  M^),  we  must,  according 

3  SO 

to  (29)  and  (30),  first  determine  .  Since  S  is  perfectly  conducting 


in  Fig.  5, 


E®  (0,  M. )  =  0  on  SSC 
— tan  —1 


where  the  subscript  tan  denotes  the  component  tangent  to  S  .  Because  the 
field  in  (31)  is  given  by  (29),  we  obtain  the  following  equation  for  J  . 


-E  (Ja,  0)  *  E  (0,  2Mn)  on  SSC 
— tan—  7  —tan  —1 


When  E  (0,  2M., )  is  evaluated  on  the  part  of  S  labeled  the  shorted 
—tan  —1 

aperture  in  Fig.  5,  the  magnetic  current  2M^  is  split  into  two  current 
sheets,  each  of  density  M^,  one  on  each  side  of  the  shorted  aperture.  As 
a  result,  E^.  (0,  2M, )  vanishes  on  the  shorted  aperture.  When  E  (0,  2M. ) 

is  evaluated  on  the  remaining  part  of  SSC,  2M^  is  taken  to  be  a  single  sheet 
of  magnetic  current  of  density  2M^  on  the  shorted  aperture. 

The  electric  current  J  in  (32)  is  approximated  by 


IN 

J8  -  I 

j=l  3  2 


where  {Ia}  are  unknown  constants  and  { }  are  the  expansion  functions  that 
appear  in  (22).  Expression  (33)  is  substituted  into  (32)  and  then  the 
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scalar  product  of  (32)  with  J  is  integrated  around  the  contour  of  S  in 


N,  we  obtain  N  simultaneous  equations 


the  xy  plane.  Letting  i  =  1,2 . 


which  can  be  written  in  matrix  form  as 


-*a 

2  I  *  v 


where  Z  is  the  PN  matrix  that  appears  in  (23)  and  each  of  the  quantities 

I3  and  V3  is  a  column  vector  of  N  elements.  The  jth  element  of  fa  is  the 

a  -^a  a 

constant  I  that  appears  in  (33).  The  ith  element  of  V  is  called  and 


is  given  by 

CN+1 

Va  =  J.  •  E(0,  2M  )dt 

x  j  1  x 

0 

where  t  and  t^+^  are  the  same  as  in  (25). 

After  the  elements  of  Z  and  V3  have  been  obtained,  (34)  can  be 

,  rj  a 

solved  for  I  .  The  electric  current  J_  will  then  be  given  by  (33)  in 

which  I.  is  the  jth  element  of  1  .  Substitution  of  (33)  into  (29)  and 
J 

(30)  gives 


L'* 

Ea(0,  tO  -  E(0,  2M  )  +  l  I3  E(J . ,  0)  in  region  a  (36) 

j-1  3  3 

and 

N 

Ha(0,  M  )  =  H(0,  2M )  +  l  Ia  H(J. ,  0)  in  region  a  (37) 

1  j-1  3  3 

The  right-hand  sides  of  (36)  and  (37)  are  the  fields  everywhere  in  Fig.  5, 
whereas  the  left-hand  sides  of  (36)  and  (37)  are  defined  only  in  region  a 
so  that  (36)  and  (37)  can  be  valid  only  in  region  a.  Substituting  (37) 
into  (11),  we  obtain 


Mj.  •  H(0,  2Mx)dy  -  l  Ij  ^  •  H~ ( J ^ ,  0)dy 

^  ^  _ T-7 


where  the  superscript  denotes  magnetic  field  evaluation  at  x-0  ,  that  is 


immediately  to  the  left  of  the  shorting  strip. 

In  (16)  and  (17),  (Eb(0,  M^),  Hb(0,  M^))  is  the  electromagnetic 
field  in  region  b  due  to  the  magnetic  current  sheet  placed  on  the 

right-hand  side  of  the  shorted  aperture  and  radiating  in  the  presence  of 

SC  sc 

S  .  Since  the  perfectly  conducting  surface  S  completely  isolates 

region  a  from  region  b,  the  magnetic  current  sheet  can  be  placed  on 
the  left-hand  side  of  the  shorting  strip  without  affecting  the  electro¬ 
magnetic  field  (]Jb(0,M^),  Hb(0,  M^))  in  region  b.  Introduction  of  this 
current  sheet  gives  the  situation  of  Fig.  5.  Since  the  electromagnetic 
field  everywhere  in  Fig.  5  is  given  by  the  right-hand  sides  of  (36)  and 
(37),  it  is  evident  that 

N 

E  (0,  MJ  =  E(0,  2M,)  +  T  I3  E(J . ,  0)  in  region  b  (39) 

-  -1  —  —L  3  3 

and 

N 

H  (0,  )  =  H(0,  2NL )  +  J  Ia  H(J. ,  0)  in  region  b  (40) 

Substituting  (40)  into  (12),  we  obtain 


w 


w 


Yb  =  - 


-w 


Mj  •  H(0,  2M1)dy  -  j  I®  J  Kj.  •  S+(Jj ,  0)dy  (41) 

3  ^  _ T.T 


-w 


where  the  superscript  "+"  denotes  magnetic  field  evaluation  at  x=0  , 
that  is,  immediately  to  the  right  of  the  shorting  strip. 

The  electromagnetic  field  (E,  H)  of  Fig.  1  is  now  given  by  (3)  in 


which  j£a  is  obtained  by  substituting  (26)  and  (36)  into  (14),  H  is  ob¬ 


tained  by  substituting  (27)  and  (37)  into  (15),  _E  is  obtained  by  substi¬ 


tuting  (39)  into  (16),  and  H°  is  obtained  by  substituting  (40)  into  (17). 


C /„inc 


(Einc  +  1  I^ncE(J.,  0))  +  V(E(0,  2M  )  +  1  iaE(J.,  0)),  region  a 

i=l  3  3  i=l  3  3 


-V(E(0,  2M. )  +  l  I®  E(J  ,  0))  ,  region  b 

j"l  J  3 

.  N  N 

(H  +  l  IincH(J.,  0))  +  V(H(0,  2M  )  +  l  laH(J.,  0)),  region 
j=l  3  3  j=l  3  3 


H  =  ( 


«  L'i 

/  -V (H(0,  2M )  +  l  I®  H(J  ,  0)) 

L  j-i  3  3 


,  region  b 


2.nc  ->inc 

In  (42)  and  (43),  1^  is  the  jth  element  of  the  column  vector  I  that 

satisfies  (23)  in  which  the  ijth  element  of  Z  is  given  by  (24)  and  the  ith 

element  of  Vinc  is  given  by  (25).  Furthermore,  I®  is  the  jth  element  of 

the  column  vector  I  that  satisfies  (34)  in  which  the  ijth  element  of  Z  is 

given  by  (24)  and  the  ith  element  of  V  is  given  by  (35).  Moreover,  V  is 

3.  b 

the  complex  number  that  satisfies  (10)  in  which  Y  ,  Y  ,  and  I  are  given  by 
(38),  (41),  and  (28),  respectively. 


VI.  EVALUATION  OF  THE  ELECTRIC  CURRENT  COEFFICIENTS  IlnC  AND  if  FOR  TM 


PLANE  WAVE  EXCITATION 

In  this  section,  the  coefficients  {l*nC^  and  rI®}  that  appear  in 

(42)  and  (43)  are  evaluated  for  the  case  in  which  the  incident  electro- 
l np  inc 

magnetic  field  (E  ,  H  )  is  the  TM  plane  wave  given  by 

.  /  Ainc  ,  .  -inc . 

^mc  _  ^  gjk(x  cos  o>  +  y  sm  ^ 


,  .inc  1  ,  .  ,inc  ,inc. 

H  =  —  (— u  sin  <t>  +  u  cos  )e 

—  U  — x  — y 


,  .inc  .  .  ,inc. 

jk(x  cos  cp  +  y  sm  <p  ) 


where  u  ,  u  ,  and  u  are  the  unit  vectors  in  the  x,y,  and  z  directions, 
—x  — y  — z 


respectively.  Moreover,  n  and  k  are,  respectively,  the  intrinsic  impedance 

and  wave  number  of  the  homogeneous  medium  with  permeability  y  and  permit- 
.  .  inc 

tivxty  E.  Furthermore,  <6  is  the  azimuthal  angle  of  the  direction  from 

•  inC 

which  the  plane  wave  comes.  Measured  in  the  xy  plane,  5  is  zero  on  the 

positive  x  axis  and  increases  as  the  y  axis  is  approached, 
inc  inc 

When  (j[  ,  ft  )  is  the  TM  field  given  by  (44)  and  (45),  we  let  the 

electric  current  expansion  function  J.  be  the  z  directed  pulse  of  unit 
amplitude  on  the  jth  segment  of  the  contour  of  S  in  the  xy  plane. 


-j(t)  = 


u  ,  t.  <  t  <  t... 
-z  j  -  j+1 


0  ,  otherwise 


,  j=l,2, . . .  ,N 


In  (46),  tj  is  the  value  of  t  at  the  point  .  The  previous  statement  is 
valid  for  j=l,2,...,N.  As  in  (24),  t^+^  is  the  total  length  of  the  contour 

SC 

of  S  in  the  xy  plane.  For  j=l , 2 , . . . ,N-1 ,  the  jth  segment  of  the  contour 

SC 

of  S  in  the  xy  plane  is  the  straight  line  segment  that  runs  from  the  point 

SC 

Pj  to  the  point  P_.+^  in  Fig-  3.  The  Nth  segment  of  the  contour  of  S  in 
the  xy  plane  is  the  straight  line  segment  that  runs  from  the  point  P^  to 
the  point  P^  in  Fig.  3. 

When  (E^nc,  H^nC)  is  the  TM  field  given  by  (44)  and  (45),  we  let  the 
magnetic  current  expansion  function  M^  be  the  function  of  y  given  by 


Mx(y)  * 


Uy  yi  -  (y/w)  ,  -w  <  y  <  w 


I  y  I  >  w 


To  solve  (23)  for  I  ,  we  must  first  evaluate  the  elements  of  Z 
and  Vinc.  The  electric  field  E (J^ ,  0)  that  appears  in  expression  (24) 
for  Z^'  is  given  by  [8,  eq.  (3-5)] 


E(J.,  0)  *  -  "2“  J.(t’)  (k^(x-x')“  +  (y-y’)  )dt'  (48) 

J  L°  J  J 

0 

(2) 

where  H  is  the  Hankel  function  of  the  second  kind  of  order  zero.  Further- 
o 

more,  (x,y)  are  the  rectangular  coordinates  of  the  point  at  which  IS (J^ ,  0)  is 
evaluated,  and  (x',y')  are  the  rectangular  coordinates  of  the  point  on  the 
contour  of  S  at  which  the  arc  length  is  t'.  In  view  of  (46),  substitution 
of  (48)  into  (24)  gives 

Ci+1  S+l 

Z  =  ~  dt  dt*  H<2)  (k^x-x')2  +  (y-y')2)  (49) 


where  (x,y)  are  the  rectangular  coordinates  of  the  point  on  the  contour  of 

SC 

S  at  which  the  arc  length  is  t. 

If  new  variables  of  integration  u  and  u'  are  defined  by 

u  =  k  (t  -  0.5  (^  +  t^+1))  (50) 


then  (49)  becomes 


In  (52), 


u’=  k  (t'  -  0.5  (tj  +  tj+1)) 


7  =  — — 

ij  4k 


°>5Yi 

0.5y 

du 

-0.5Yi 

-0.5y\ 

du'  H(2)(a. .) 
o  ij 


y±  -  *  /<* i+1  -  *i>2  +  <yi+i  -  yi>2 


where  (x^,  y.)  are  the  (x,y)  coordinates  of  the  point  on  the  contour  of 
SSC  at  which  the  arc  length  is  t..  Furthermore, 


(kx .  -  kx .  +  u  cos  4>.  -u'cos  <J>.)  +  (ky.  -  ky  +  u  sin  <p  -  u  sin  4>.) 


where  x^  and  y^  are,  respectively,  the  x  and  y  coordinates  of  the  midpoint 
of  the  ith  segment  of  the  contour  of  SSC,  and  <J>  is  the  angle  that  the  ith 
segment  makes  with  the  positive  x  axis.  More  precisely, 


x.  =  0.5(x1+1  +  xA) 


yi  =  °*5(yi+l  +yi> 


k(xi+i  -  x.) 

cos  =  - 

1  Yi 


sin  4> .  = 
1 


k(Y+i  -  V 


If  i  4  j ,  the  double  integral  in  (52)  is  evaluated  by  the  method 
of  two-dimensional  Gaussian  quadrature.  If  i  *  j,  (52)  reduces  to 


0.5y.  0.5Y. 

t  i  r  i 


Z  =  — J- 
ii  4k 


du’  H^2)(|u-u'|) 


-0 . 5y .  -0 . 5y . 

i  i 


The  bracketed  term  on  the  right-hand  side  of  the  identity 


H<2)(x)  -  1  -  &  In  £p)  +  [H^2)(x)  -  1  +  &  In  (^)]  (6 

O  7T  Z  O  7T  Z 

approaches  zero  as  x  approaches  zero  [9,  p.  462].  Here,  In  y  is  Euler's 


constant.  Substitution  of  (60)  into  (59)  gives 


l  .  r- 


0.5yi 

du  f  du'  [H(2)(|u-u' I) 


!  +  !n  (112=aJ-)] 


25 


0  5Yi 

0.5Y 

du 

j 

f 

J 

-0.5Y. 

l 

-0.5y 

Using  [10,  Formulas  610.  and  610.1.]  to  evaluate  the  double  integral  in 
(62),  we  obtain 

.  2 
2  JYi 

C  *  Yj  +  -^  (3  -  2  In  (~))  ( 


Substitution  of  (63)  into  (61)  gives 

n  )  2  jYi  YYi 

zii  =  ijYi+-ir  <3-21n  (_r)) 

Q.SYi  0.5y.  _  (  1 

+  |  du  du'  [H^2)  ( |u-u'  | )  -  1  In  (yJ.4T.u'.L)]  1  (64) 

-0.5Yi  -0.5Y.  J 

Since  its  integrand  is  bounded,  the  iterated  integral  in  (64)  can  be 
accurately  evaluated  by  the  method  of  two-dimensijnal  Gaussian  quadrature. 
Substituting  (44)  and  (46)  into  (25),  we  obtain 

Ci+1  .  .  . 

,  !  ,  .me  .  .  ,inc 

„inc  jk(x  cos  <p  +  y  sin  d> 

V.  =  eJ  dt  (65) 


where  (x,y)  are  the  rectangular  coordinates  of  the  point  on  the  contour 
SC 

of  S  at  which  the  arc  length  is  t.  If  the  variable  of  integration  in 
(65)  is  changed  from  t  to  u  of  (50),  then  (65)  becomes 


,  +  ,inc  ,  +  ,  ,inc.  0.5y.  .  ,,  Ainc. 

^  jk(x^  cos  +  y^  sin  4>  )  r  i  ju  cos  ) 


~0 . 5y . 


The  integral  in  (66)  being  tractable,  (66)  reduces  to 


w  Y,  jk(x  cos  cT  +  y,  sin  <T  )  sin  (0.5y,  cos  (<{>.-<{>)) 

=  -rr  e  1  1  ( - — - 1 - ) 

O.Sy^cosC^^-^^"110) 


Now,  the  quantities  Z  and  V  that  appear  in  (23)  are  known,  the  diagonal 
elements  of  Z  being  given  by  (64),  the  remaining  elements  of  Z  by  (52),  and 
the  elements  of  V^nc  by  (67).  Next,  (23)  is  solved  for  I^nc.  The  coef¬ 
ficients  {linC}  will  then  be  the  elements  of  Ilnc. 

J 

To  solve  (34)  for  I  ,  we  must  first  evaluate  the  elements  of  V  . 

In  the  E1(0,  2M^)  that  appears  in  (35),  2M^  represents  the  two  magnetic  cur¬ 
rents,  each  of  density  M^,  that  straddle  the  y  axis  in  Fig.  5.  Therefore, 
_E(0,  2M^)  vanishes  on  the  y  axis.  Because  is  on  the  y  axis  and  because 
_E(0,  2M^)  vanishes  there,  it  is  evident  from  (35)  that 

V*  -  0  (68) 

The  next  two  paragraphs  are  devoted  to  evaluating  for  i  ^  N. 

The  electric  field  E:(0,  2M^)  that  appears  in  expression  (35)  for 
V*  is  given  by  [9,  eqs.  (3-4)  and  (5-96)] 


E(0,  2^)  =  ±  V  x  ^(y’)  H*  J(k  /  x  +  (y-y’)  )  dy’ 


where  (x,y)  are  the  rectangular  coordinates  of  the  point  at  which  _E(0,  2M^) 
is  evaluated.  Furthermore,  "V*"  is  the  curl  operator  with  respect  to  the 
coordinates  (x,y).  If  i  ?  N,  then  _J  is  remote  from  the  magnetic  current 
sheet  2M^  so  that  the  E(0,  2M^)  that  appears  in  (35)  is  not  evaluated  on 
this  magnetic  current  sheet.  In  this  case,  the  curl  operator  in  (69)  can 


be  taken  under  the  integral  sign.  Taking  the  curl  operator  under  the 
integral  sign,  substituting  (47)  for  in  (69),  and  using  [9,  eqs.  (A-14) 
and  (D-15)],  we  obtain 


(y'/w)  H;  (kj/x  +  (y-y  )  ) 


-w  y  x  +  (y-y') 

Substitution  of  (70)  and  (46)  into  (35)  gives 


"r+1  wr  x\l  1  -  (y'/w)2  Hf2)(kv/x2  +  (y-y')2) 


t .  -w 
l 


+  (y-y') 


where  (x,y)  are  the  rectangular  coordinates  of  the  point  on  the  contour  of 

SC  i 

S  at  which  the  arc  length  is  t.  If  new  variables  of  integration  u  and  u 
are  defined  by 


then  (71)  becomes 


u  =  k(t  -  0.5  (ti  +  ti+1)) 


u  *  =  ky  ' 


,  #fSY*n  +  1  n  ,r  nl 

i  2k  j  i  i  J 

-°*5Yi  -kw  (74) 


where  Y^,  cos  and  a^N  are  given  by  (53),  (55),  (57),  and  (54),  respec 

SC 

tively.  Because  of  the  orientation  of  the  Nth  segment  of  the  contour  of  S  , 
expression  (54)  for  a...  reduces  to 


(kx^  +  u  cos  $^)  +  (ky^  +  u  sin  <p^  -  u  )  (75) 

The  double  integral  in  (74)  is  evaluated  by  the  method  of  two-dimensional 
Gaussian  quadrature.  Now,  the  quantities  Z  and  V  that  appear  in  (34)  are 
known,  the  diagonal  elements  of  Z  being  given  by  (64),  the  remaining  elements 


of  Z  by  (52),  V^  by  (68),  and  the  remaining  elements  of  V  by  (74).  Ne 
(34)  is  solved  for  Ia.  The  coefficients  {if}  will  then  be  the  elements 
of  la. 


VII.  EVALUATION  OF  THE  MAGNETIC  CURRENT  COEFFICIENT  V  FOR  TM  PLANE 


WAVE  EXCITATION 

In  this  section,  the  coefficient  V  that  appears  in  (42)  and  (43) 

is  evaluated  for  the  case  in  which  the  incident  electromagnetic  field 
inc  inc 

(E  ,  H  )  is  the  TM  plane  wave  given  by  (44)  and  (45),  the  electric  cur¬ 
rent  expansion  functions  {Jj }  are  given  by  (46),  and  the  magnetic  current 
expansion  function  is  given  by  (47).  To  solve  (10)  for  V,  we  must  first 
evaluate  Ya,  Yb,  and  I.  Here,  Ya  is  given  by  (38),  Yb  by  (41),  and  I  by 

g 

(28).  The  next  eight  paragraphs  are  devoted  to  evaluating  Y  . 

£ 

Expression  (38)  for  Y  is  recast  as 


where 


Ya  -  Yhs  -  l  I*  CT 

j-1  J  J 


•  H(0,  2M,  )dy 


Mi  *  h  (2-j* 


The  superscript  "hs"  in  (77)  denotes  half  space.  As  defined  by  (77),  Y  is 

g 

what  Y  would  be  if  region  a  were  the  half-space  region  for  which  x  <  0.  The 
truth  of  this  statement  is  evident  because,  if  region  a  were  the  half-space 

SC 

region  for  which  x  <  0,  then  S  would  be  the  yz  plane  at  x*0,  and  the  split 


magnetic  current  sheet  2M^  defined  in  the  sentence  that  contains  (29)  would 
induce  no  electric  current  on  this  plane  so  that  the  summation  with  respect 


to  j  in  (76)  would  vanish. 

The  H(0,  2M^)  that  appears  in  (77)  is  obtained  by  replacing  _ES  by 
H(0,  211^),  J  by  v  by  e,  and  £  by  u  in  [8,  eqs.  (3-52)  to  (3-54)].  Thus 


where 


H(0,  2M  )  .-joiF  -  V<f> 


jr  J  Mx (y ' )  H^2)(k]y-y'|)dy’ 


i  7^Vi2h(2) 


*  -^T-Hr^iy-y’Ddy' 


Here,  u)  is  the  angular  frequency.  It  has  been  assumed  that  M^(y')  exists 
only  at  x  =  0  and  has  only  a  y  component  M^(y').  The  (x,y)  coordinates 
at  which  F  and  $  are  evaluated  in  (80)  and  (81)  are  (0,  y).  Substitution 
of  (79)  into  (77)  gives 


h  s 

Y  =  jco  (y)  •  Fdy  +  tt^ty)  dy  dy 


Integrating  by  parts  the  second  integral  on  the  right-hand  side  of  (82)  and 
anticipating  that  M^(±w)  will  vanish,  we  obtain 

hs  T  7  d 

Y  *  j  to  M.  (y)  •  Fdy  -  - -r- -  $  dy  (83 

J  -i  —  J  dy 

-w  -w 

If  (80),  (81),  and  (47)  are  substituted  into  (83)  and  if  new  variables  of 
integration  u  and  u'  are  defined  by 


then  (83)  becomes 


Tv  rYKV^vjrVFTjJTV^  ^  VWI  M.Vflnw  -n 


M2  rfi,  ’  «  — 


-  (u')  du 


From  [10,  Formula  352.01.1,  we  obtain 


2  h  2  ,  it 
u  yl  -  u  du  -  g- 


Equations  (93)  -  (95)  are  used  to  reduce  (91)  to 


_  tt  .  iir  ,1  ,  ,ykws. 

h  •  T  +  2  <4  -  ln  » 


It  is  evident  from  [11,  Eq.  (A6b)]  that 


ln I u-u '  du '  * 


fl  -  (u'V 


It  is  evident  from  [10,  Formula  322.01.]  that 


.  w 

du  =  y 


-1  ,/ 1  ~ 


Equations  (97)  and  (98)  are  used  to  reduce  (92)  to 


I*  -  Jir  W 

l.  g  ^ 

Now,  the  Y  that  appears  in  expression  (76)  for  Y  is  given  by 

(86)  in  which  C.  and  C,  are  given  by  (89)  and  (90),  respectively.  In  (89) 
A  v 

and  (90),  I  and  I$  are  given  by  (96)  and  (99),  respectively.  The  double 
integrals  that  appear  explicitly  in  (89)  and  (90)  are  evaluated  by  the 
method  of  two-dimensional  Gaussian  quadrature. 

The  C~  that  appears  in  expression  (76)  for  Ya  is  given  by  (78).  In 
(78),  the  y  component  of  H  (J^,  0)  is  called  R^OJ^,  0).  Since  Is  u^ 


on  the  shorting  strip,  it  follows  that,  at  x*0~, 


V4- 


0) 


r-o.5, 

0  ,  i  y  |  >  w 


Substitution  of  (47)  and  (100)  into  (78)  gives 


(100) 


cs  -  -  T  I  i/1  -  <£>2 


(101) 


Expressing  the  integral  in  (101)  as  the  product  of  w  with  the  integral  in 
(94) ,  we  obtain 


TTW 

4 


(102) 


Striving  to  evaluate  C^.  for  j  4  N,  we  let  the  y  component  of  H(J^ ,  0) 
be  Hy(J^,  0)  and  express  it  as  [8,  eqs.  (3-26)  and  (3-32)] 


N+l 

H  (J.,  0)  =  u  •  7*  ~  f  J.(t')  H(2)  (k  (x— x 1 ) 2  +  (y-v ' ) 2)dt ' 

y  -J  -y  4j  j  -3  o  v 


(103) 


where  (x,y)  are  the  rectangular  coordinates  of  the  point  at  which  Hy(J^,  0) 
is  evaluated,  and  (x',  y')  are  the  rectangular  coordinates  of  the  point  on 
the  contour  of  S  at  which  the  arc  length  is  t'.  In  (103),  "V  x"  is  the 
curl  operator  with  respect  to  the  coordinates  (x,y).  If  (x',  y')  never 
coincides  with  (x,y),  then  the  curl  operator  can  be  taken  under  the  integral 
sign  in  (103).  Taking  the  curl  operator  under  the  integral  sign,  substituting 
(46)  into  (103),  and  using  [9,  eqs.  (A-14)  and  (D-15)],  we  obtain 


(J  0)  .  _jk  (x-x')  (k  ^(x-x'T  +  (y-y ')  ) 


2  2 

(x-x')  +  (y-y*) 


dt’  (104) 


Substitution  of  (47)  and  (104)  into  (78)  gives 

-  ik  1  ! - I  V1  dt'x'  H{2)<k/(x')2  +  (y-y')2) 

c.  .  f  jdy  y  l-(y/w;  —^.1^=7. -  >  i*  » 

-w  t  +  (y-y')  (105) 

If  new  variables  of  integration  u  and  u'  are  defined  by 


u  =  ky 


(106) 


u'=  k(t’  -  0.5(t  +  t.+1)) 


then  (105)  becomes 


c  —  J- 
j  4k 


(107) 


r — 7772  °rSYj  *,,(l“t  +  u'  cos  V  hi2)(cSi1)  . 


-0.5y. 


(108) 


where  y  ,  x^  ,  cos  <J>^ ,  and  are  given  by  (53),  (55),  (57),  and  (54), 
respectively.  The  double  integral  on  the  right-hand  side  of  (108)  is  the 
double  integral  on  the  right-hand  side  of  (74)  with  i  replaced  by  j,  with  u 
and  u'  interchanged,  and  with  the  order  of  integration  interchanged. 
Therefore, 


j-1,2 . N-l 


(109) 


where  Va  is  given  by  the  right-hand  side  of  (74)  with  i  replaced  by  j. 

Now,  the  {Cj}  that  appear  in  expression  (76)  for  Ya  are  given  by 

h  s 

(109)  and  (102).  The  evaluation  of  the  Y  that  appears  in  (76)  was  de¬ 
scribed  in  the  paragraph  that  follows  (99).  The  coefficients  {la}  that 
appear  in  (76)  were  evaluated  in  Section  VI.  This  concludes  our  de- 
scription  of  the  evaluation  of  Y  . 


Expression  (41)  for  Y  is  recast  as 


,b  „hs 


=  y  -  y  I  C 
J-l  J  3 


tlS  <f 

where  Y  is  given  by  (77)  and  is  given  by 

w 

C+  =  I  M.  <  H+(J.,  0)  dy  (1 

J  J  ”  -3 

-w 

"4" 

In  (111),  the  y  component  of  H  (J.Y,  0)  is  called  H  (J.Y,  0).  Since  J.,  is 

—  — N  y  — N  — N 

uz  on  the  shorting  strip,  it  follows  that,  at  x  =  0+ 


V4- 0)  - 


0.5  ,  -w  <  y  <  w 


0  ,  i  y  |  >  w 


Now,  H*  of  (112)  is  the  negative  of  of  (100).  Therefore,  C+  of  (111) 
is  the  negative  of  CN  of  (102). 

<  •  -  (1 

If  j  N,  then  C*  of  (111)  is  equal  to  of  (109). 

cT  =  cT»  j  =  1.2.....N-1  (1 

J  j 

Now  the  t,C+}  that  appear  in  expression  (110)  for  Y*3  are  given  by 
(114)  and  (113).  The  quantities  Y^S  and  I*  on  the  right-hand  side  of 
(110)  are  the  same  as  in  (76).  This  concludes  our  description  of  the 
evaluation  of  Y*3. 

The  quantity  (Y3  +  Y*5)  appears  In  (10).  Thanks  to  (109),  (113), 
and  (114),  the  sum  of  (76)  and  (110)  reduces  to 
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Y3  +  Yb  =  2Yhs  +  '  T  Ia  Va 


j-1 


J  J 


(115) 


Expression  (28)  for  I  is  recast  as 


i  =  clnc  +  l  IincC“ 

j-1  3  J 


(116) 


where  C_.  is  given  by  (78),  and 


w 


Clnc  =  |  M,  •  H1 

-w 


Mx  •  Hlncdy 


(117) 


Substituting  (45)  and  (47)  into  (117)  and  then  introducing  the  new  variable 


of  integration  u  defined  by  (84),  we  obtain 

1 


„inc  w  .inc 

C  =  —  cos  d> 


^1  -  u2  cos  (kw  u  sin  $inc)du 


(118) 


-1 


The  integral  in  (118)  is  evaluated  by  the  method  of  Gaussian  quadrature. 

Now,  the  quantity  I  that  appears  in  (10)  is  given  by  (116)  in  which 


ClnC  is  given  by  (118),  the  tCj}  are  given  by  (109)  and  (102),  and  the 


inc  a  b 

(Ij  }  were  evaluated  in  Section  VI.  The  quantity  (Y  +  Y  )  that  appears 


hs  . 


in  (10)  is  given  by  (115)  in  which  Y  is  evaluated  as  described  in  the 


paragraph  that  follows  (99),  the  {i^}  were  evaluated  in  Section  VI,  and 


a  a  b 

V  is  given  by  (74)  with  i  replaced  by  j.  After  I  and  (Y  +  Y  )  have  been 


evaluated,  (10)  can  be  solved  for  V. 


VIII.  NUMERICAL  RESULTS  FOR  A  SLOTTED  TM  CIRCULAR  CYLINDRICAL  SURFACE 


In  this  section,  numerical  results  are  presented  for  the  magnetic 
current  coefficient  V  for  the  slotted  TM  cylindrical  surface  shown  in 


Fig.  6  with 


ka  =  t/2 


4>  =5* 
o 


(119) 


xnc 

In  (44),  $  is  set  equal  to  180°  so  that  the  incident  electric  field 

„inc  .  .  , 

_E  is  given  by 


_inc  -jkx 

E  =  u  e 
—  — z 


(120) 


The  coefficient  V  was  calculated  by  three  different  methods,  the  method 
of  solution  with  pseudo-image,  the  method  of  solution  without  pseudo¬ 
image,  and  the  Fourier  series  method  of  solution.  The  method  of  solution 
with  pseudo-image  is  the  method  of  solution  described  in  Sections  III  to 
VII.  This  method  was  applied  using  10-point  Gaussian  quadrature  whenever 
Gaussian  quadrature  was  called  for,  dividing  the  circular  arc  of  Fig.  6 
into  35  small  arcs  each  of  angular  extent  10°,  and  replacing  each  of  these 
arcs  by  the  straight  line  segment  joining  its  end  points.  The  method  of 
solution  without  pseudo-image  is  described  in  Appendix  A,  and  the  Fourier 
series  method  of  solution  is  described  in  Appendix  B.  All  three  methods 
of  solution  were  implemented  by  means  of  computer  programs  which  will  be 
described  and  listed  in  a  forthcoming  report. 

The  method  of  solution  without  pseudo-image  differs  from  the  method 
of  solution  described  in  Sections  III  to  VII  in  that  the  pseudo-image  of 
Section  V  is  eliminated.  That  is,  the  short-circuit  fields  due  to  the 


magnetic  current  expansion  function  placed  on  the  left-hand  side  of  the 


The  slotted  TM  cylindrical  surface  as  viewed  in  the  xy  plane.  It  is 

an  infinitesimally  thin  perfectly  conducting  circular  cylindrical 

shell  of  radius  a  centered  at  (x-a  cos  4>  ,  y*0)  with  a  slot  of 

o 

angular  width  2 <J>  .  The  edges  of  the  shell  are  at  (x*0,  y*±a  sin  4>q) 
ka  *  tt/2  and  <J>  ■  5a. 
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shorted  aperture  are  obtained  by  solving  the  scattering  problem  of  Fig.  k 
rather  than  that  of  Fig.  5.  Similarly,  the  short-circuit  fields  due  to 
placed  on  the  right-hand  side  of  the  shorted  aperture  are  obtained  by  solv¬ 
ing  the  scattering  problem  not  of  Fig.  5  but  of  the  situation  that  would 
exist  in  Fig.  4  if  were  placed  on  the  right-hand  side  of  the  shorted 
aperture  instead  of  on  the  left-hand  side  of  it. 

In  the  Fourier  series  method  of  solution,  the  aperture  electric  field 
rather  than  the  magnetic  current  M  is  expanded.  Furthermore,  JE  is  ex¬ 
panded  not  on  the  straight  line  connects  the  end  points  (x  =  0,  y  =  t  a  sin  4>q) 
of  the  circular  arc  in  Fig.  6,  but  on  the  circular  arc  of  angular  extent  2 <J>q 
that  connects  these  end  points.  The  combination  of  this  circular  arc  of  angu¬ 


lar  extent  2$  and  the  existing  circular  arc  in  Fig.  6  would  form  a  complete 
o 


circle.  In  the  Fourier  series  method  of  solution. 


E(<t>)  =  -  Vu  \Jl  -  (^-~-^8--)2  ,  -<D  1 


d>  -  180°  <  <j> 
o  —  —  o 


(121) 


where  V  is  an  unknown  complex  constant  to  be  evaluated  and  <f>  is  the  angle 
shown  in  Fig.  6.  If  the  y  axis  passed  through  the  center  of  the  circle 
whose  arc  appears  in  Fig.  6,  then  4  would  be  the  usual  cylindrical  coordi¬ 
nate  angle  measured  in  the  counterclockwise  direction  from  the  positive 
x  axis.  In  (121),  E(4>)  is  the  electric  field  on  the  circular  arc  of  extent 
24  needed  to  complete  the  circle  in  Fig.  6.  Therefore,  the  (x,y)  coordi¬ 


nates  of  the  location  of  E(4)  ate  given  by 


x  =  a(cos  4  +  cos  4) 
o 


(122) 


y  *  a  sin  | 


To  compare  the  method  of  solution  with  pseudo-image  to  the  Fourier 
series  method  of  solution,  we  must  interpret  M  of  (9)  in  terms  of  the 
electric  field  £  in  the  slot.  Combining  (1)  with  (9)  and  assuming  that  E 
has  no  component  in  the  x  direction,  we  obtain 

E  =  -  V  u  *  M-  (123) 

—  — x  — 1 

Substitution  of  (47)  into  (123)  gives 

(y )  =  -  Vuz/l  -  (y/w)^  ,  -w  _<  y  £  w  (124) 

No  confusion  should  arise  from  the  fact  that  the  argument  of  E  of  (124) 

is  a  rectangular  coordinate  whereas  that  of  IS  of  (121)  is  an  angle.  The 

_E  of  (124)  is  the  electric  field  at  x  =  0  in  the  slot,  but  the  JS  of  (121) 

is  the  electric  field  on  the  circular  arc  of  angular  extent  2^  needed  to 

complete  the  circle  in  Fig.  6.  The  (x,y)  coordinates  on  this  circular 

arc  are  given  in  terms  of  (j>  by  (122) . 

Calculated  values  of  V  are  presented  in  Table  1.  In  Table  1,  the 

phases  of  V  are  nearly  the  same,  differing  by  0.24°  at  most.  The  phase 

lead  of  0.21°  for  the  V  obtained  by  the  Fourier  series  method  of  solution 

with  respect  to  the  V  obtained  by  the  method  of  solution  with  pseudo-image 

inc 

is  reasonable  because  it  is  shown  as  follows  that  the  phase  of  12  at  the 

inc 

location  of  E!(180o)  of  (121)  is  0.342°  and  that  the  phase  of  E  at  the 
location  of  E(0)  of  (124)  is  0°.  According  to  the  first  of  equations  (122) 
and  the  second  of  equations  (119),  the  x  coordinate  of  the  location  of 
E(180°)  of  (121)  is  given  by 


x 


-  0.0038a 


(125) 


inc 

Thanks  to  the  first  of  equations  (119),  the  phase  of  jS  of  (120)  is 
0.342°  when  x  is  given  by  (125).  Therefore,  the  phase  of  E*nc  at  the 
location  of  £(180°)  of  (121)  is  0.342°.  The  x  coordinate  of  the  loca¬ 
tion  of  £(0)  of  (124)  is  zero.  At  x  =  0,  the  phase  of  E^nc  of  (120)  is  0°. 

inc 

Therefore,  the  phase  of  E  at  the  location  of  E(0)  of  (124)  is  0°. 

Table  1.  The  coefficient  V  for  the  slotted  TM  circular  cylindrical 
surface  shown  in  Fig.  6. 


Method  of  solution 

Real 

part  of  V 

Imaginary 
part  of  V 

Magnitude 
of  V 

Phase  of  V 
in  degrees 

With  pseudo-image 

-0.0347 

-0.1483 

0.1523 

-103.17 

Without  pseudo-image 

-0.0399 

-0.1701 

0.1747 

-103.20 

Fourier  series 

-0.0344 

-0.1495 

0.1534 

-102.96 

In  Table  1,  the  magnitude  of  the  V  obtained  by  the  method  of  solu¬ 
tion  with  pseudo-image  is  a  trifle  smaller  than  that  of  the  V  obtained 
by  the  Fourier  series  method  of  solution.  This  is  expected  because  the 
location  of  _E(0)  of  (124)  is  slightly  inside  the  circular  arc  in  Fig.  6 
as  compared  to  the  location  of  _E(180°)  of  (121).  The  electric  field 
should  decrease  as  the  interior  of  an  aperture  perforated  cavity  is 
approached. 

In  Table  1,  the  complex  number  V  obtained  by  the  method  of  solu¬ 
tion  with  pseudo-image  agrees  well  with  that  obtained  by  the  Fourier 
series  method  of  solution.  However,  the  magnitude  of  the  complex  number 
V  obtained  by  the  method  of  solution  without  pseudo-image  is  more  than 
10%  too  large.  In  this  case,  the  complex  number  V  obtained  by  the  method 
of  solution  with  pseudo-image  is  evidently  more  accurate  than  that  obtained 
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by  the  method  of  solution  without  pseudo-image.  We  believe  that*  for 
scattering  problems  involving  small  apertures,  the  method  of  solution 
with  pseudo-image  gives  results  that  are  generally  more  accurate  than 
those  obtained  by  the  method  of  solution  without  pseudo-image.  More 
computations  are  presently  being  done  in  order  to  verify  this  belief. 
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APPENDIX  A 


THE  METHOD  OF  SOLUTION  WITHOUT  PSEUDO-IMAGE 


Calculation  of  the  magnetic  current  coefficient  V  by  means  of  the 

method  of  solution  without  pseudo-image  is  described  in  this  appendix. 

3  b 

Equation  (10)  can  be  solved  for  V  if  the  quantities  Y  ,  Y  ,  and  I  are 
known.  In  the  method  of  solution  without  pseudo-image,  the  short-circuit 
magnetic  field  H  (0,  M^)  in  expression  (11)  for  Y  is  obtained  by  solving 
the  scattering  problem  of  Fig.  4  rather  than  that  of  Fig.  5.  Moreover, 
the  short-circuit  magnetic  field  H^ (0,  M^)  in  expression  (12)  for  Y^  is 
obtained  by  solving  the  scattering  problem  not  of  Fig.  5  but  of  the 
situation  that  would  exist  in  Fig.  4  if  M^  were  placed  on  the  right-hand 
side  of  the  shorted  aperture  instead  of  on  the  left-hand  side  of  it. 
Furthermore,  the  quantity  I  in  (10)  is  the  same  as  obtained  by  the  method 
of  solution  with  p seudo- image . 

Considering  the  situation  of  Fig.  4  rather  than  that  of  Fig.  5,  we 
write,  instead  of  (29)  and  (30), 


Ea(0,  Mx)  =  E(0,  M~)  +  E(JaW,  0)  in  region  a 


Ha(0,  Mx)  »  H(0 ,  M~)  +  H(JaW,  0)  in  region  a  (A- 

where  Jaw  is  the  electric  current  induced  on  SSC  in  Fig.  4.  The  "w"  in 
the  superscript  "aw"  means  without  pseudo-image.  The  superscript  "-"  on 
M^  on  the  right-hand  sides  of  (A-l)  and  (A-2)  indicates  that  M^  is  placed 

3W 

on  the  left-hand  side  of  the  shorted  aperture.  Analogous  to  (33),  J  is 


approximated  by 


so  that  (34)  is  replaced  by 


-*-aw 
Z  I 


vaw 


(A-4) 


where  the  jth  element  of  I  W  is  the  constant  1^  that  appears  in  (A-3). 

-►aw 

The  ith  element  of  V  is  given  by 


Sl+l 

J±  •  E(0,  M")dt 

0 


(A-5) 


Substitution  of  (A-3)  into  (A-2)  gives 


Ha(0,  Mn)  -  H(0,  mT)  +  1  if*  E(J.»  0)  in  region  a  (A-6) 

1  J*1  2  J 

a  aw 

Substituting  (A-6)  into  (11)  with  Y  replaced  by  Y  where  the  extra 

superscript  "w"  means  without  pseudo- image,  we  obtain 

.  I  Yhs  -  f  I "  cT  (A- 7 

2  j-1  7  J 

Vis  — 

where  Y  is  given  by  (77)  and  by  (78). 

Consider  the  situation  that  would  exist  if  was  on  the  right-hand 

side  of  the  shorted  aperture  in  Fig.  4  instead  of  on  the  left-hand  side  of 

it.  In  this  situation,  the  electromagnetic  field  (Eb(0,  ML^),  Hb(0,  M^))  is 

given  by 


Eb(0,  l^)  -  E(0,  l£)  +  E(JbW,  0) 

in  region  b 

(A-8) 

Hb(0,  Mx)  -  H(0,  M*)  +  H(JbW,  0) 

in  region  b 

(A— 9 ) 

where  jbw  is  the  electric  current  induced  on  SSC  in  the  previously  de¬ 
scribed  situation.  In  the  superscript  "bw"  the  "b"  denotes  the  electro¬ 
magnetic  field  problem  in  region  b,  and  the  "w"  means  without  pseudo-image. 


The  superscript  "+"  on  on  the  right-hand  sides  of  (A-8)  and  (A-9)  indi 

cates  that  is  placed  on  the  right-hand  side  of  the  shorted  aperture. 

bw 

Analogous  to  (33) ,  J[  is  approximated  by 


_bw 


N  . 

I  TbWJ. 
j-1  J 


(A- 


Using  the  method  of  moments  to  solve  the  electromagnetic  field  problem, 
we  obtain  the  following  matrix  equation  which  is  similar  to  (34) . 


-►bw  -*-bw 

z  r  -  V 


(A- 


The  jth  element  of  Ibw  is  the  constant  Ibw  that  appears  in  (A-10) . 
-*bw 

ith  element  of  V  is  given  by 


The 


N+l 


,bw 


*  E(0,  M1)dt 


(A- 


Substitution  of  (A-10)  into  (A-9)  gives 


N  , 

H  (0,  M.)  -  H(0,  M-  )  +  l  I.”  H(J  , 

j-1  J  3 


0) 


in  region  b 


(A- 


b  bw 

Substituting  (A-13)  into  (12)  with  Y  replaced  by  Y  where  the  extra 
superscript  "w"  means  without  pseudo-image,  we  obtain 


Ybv .  1  T1»  .  2 

j-1  3  3 

where  Ybs  is  given  by  (77)  and  by  (111) . 
The  sum  of  (A-7)  and  (A-14)  is 


Yaw  +  ybw  ,  Yhs 


N 


-  £  + 

*1"  — 

Because  is  given  by  (102)  and  (109)  and  is  given  in  terms  of 


(A- 


(A- 


M 


tf  j  -r  ■ 


by  (113)  and  (114) ,  (A-15)  becomes 


Yaw  +  Ybw  =  yhs  ‘  TTW  '-aw 


N-l 


+  ™  ^  +  7  l  ar'  +  ifiv* 

]-l  ]  3  3 


4  '~N  ‘N 
The  sum  of  (A-4)  and  (A-ll)  is 


(A-H 


z[iaw  +  ibw]  -  vaw  +  vbw 


(A-l 


It  is  evident  from  (A-5),  (A-12),  and  (35)  that 


->aw  -*bw  -*a 
V  +  v  *  V 


(A-ll 


Substituting  (A-18)  into  (A-17)  and  comparing  the  result  with  (34),  we 
obtain 


+  IbW 
J 


j-1.2. 


,N 


(A-l1 


Equation  (A-19)  is  true  because  the  sum  of  the  electromagnetic  sources  in 
the  field  problems  whose  matrix  equations  are  (A-4)  and  (A-ll)  is  the  elec¬ 
tromagnetic  source  in  the  field  problem  whose  matrix  equation  is  (34). 
Subtracting  (A-ll)  from  (A-4),  we  obtain 


-►aw  -*bw,  ■►aw  **bw 
Z[I  -  I  ]  =  V  -  V 


(A-2 


The  difference  between  (A-5)  and  (A-12)  is 

t. 


vaw  -  vbw 

i  i 


N+l 

. 


J.  •  (E(0,  Ml)  -  E(0,  M1))dt 


(A-2 


Because  M  is  ^  placed  on  the  left-hand  side  of  the  shorted  aperture 
and  M+  is  placed  on  the  right-hand  side  of  the  shorted  aperture  where 
is  given  by  (47),  the  difference  field  in  (A-21)  exists  only  at  x  *  0 
between  the  magnetic  current  sheets  M.  and  where  it  is  given  by 


Substituting  (A-22)  and  (46)  into  (A-21)  and  using  (94),  we  obtain 


Substitution  of  (A-19)  into  (A-16)  gives 


(A- 


Yaw  +  ybW  =  YhS 


irw  .  aw  _  ±  y  »  v- 

4  LN  ;  2  *  J  j 


„bw 


,  N-l 

1  r  ,a  „a 


(A- 


where,  from  what  was  done  in  the  previous  paragraph,  (IaW  ~  IbW)  is  the 


Nth  element  of  the  solution  [I3W  -  Ibw]  to  (A-20)  when  the  elements  of 


[VaW  -  V  W]  are  given  by  (A-23) .  Comparison  of  (A-24)  with  (115)  gives 


yaw  +  ybw 


1  +  Yb,  +  f  (if  -  l£»> 


(A- 


Affixing  the  superscript  "w"  to  all  quantities  on  the  left-hand  side  of 

(10)  to  indicate  without  pseudo-image,  and  solving  the  resulting  equation 
w 

for  V  ,  we  obtain 

VW  =  l/(YaW  +  YbW)  (A- 


where  VW  is  the  magnetic  current  coefficient  obtained  by  the  method  of 
solution  without  pseudo-image  and  (Yaw  +  Ybw)  is  given  by  (A-25). 


APPENDIX  B 


THE  FOURIER  SERIES  METHOD  OF  SOLUTION  FOR  A  SLOTTED  TM 


CIRCULAR  CYLINDRICAL  SURFACE 


In  this  appendix,  the  Fourier  series  method  of  solution  is  used  to 

inc 

solve  the  scattering  problem  that  arises  when  the  electric  field  E  of 

(120)  is  incident  on  the  conducting  surface  in  Fig.  6.  We  require  that 

E^  be  continuous  at  p*a  where  is  the  z  component  of  the  electric  field 

and  p  is  the  distance  from  the  center  of  the  circle  whose  arc  appears  in 

Fig.  6.  Combining  (121)  with  the  fact  that  E  vanishes  on  the  circular 

z 

arc  in  Fig.  6,  we  obtain 


-  V  y  1  -  (*  "y  -0?  ,  -0o  £<f>  -  180°  <  <J>o 


,  other  values  of  <p 


(B-l) 


at  0=a.  If  the  y  axis  passed  through  the  center  of  the  circle  whose  arc 

appears  in  Fig.  6,  then  4>  would  be  the  usual  cylindrical  coordinate  angle 

measured  in  the  counterclockwise  direction  from  the  positive  x  axis.  Since 

E  is  constrained  by  (B-l),  the  4>  component  of  the  magnetic  field  can  not  be 
z 

continuous  across  the  entire  circular  arc  of  angular  extent  24>q  that  bridges 
the  gap  to  form  a  complete  circle  in  Fig.  6.  However,  the  4>  component  of  the 
magnetic  field  can  satisfy 


180°+4> 

j  °  <«*  -  V  i/1  -  c4^21)2  d* ' 0 

180-4)  ° 
o 


(B-2) 


where  H~  is  the  limit  as  p  approaches  a  from  below  of  the  4>  component  of 

4> 

the  magnetic  field  and  H*  is  the  limit  as  p  approaches  a  from  above  of  the 


I  *  *  •  ' 


4>  component  of  the  magnetic  field.  In  what  follows,  E  is  expanded  in 

z 

separate  Fourier  series  in  the  regions  for  which  0  <  a  and  c  >  a,  E  is 

z 

required  to  be  continuous  at  p*a,  and  E  is  forced  to  conform  with  (B-l) 

— z 

at  c*a.  Finally,  (B-2)  is  used  to  determine  the  unknown  complex  constant 
V  that  appears  in  (B-l). 

In  the  region  for  which  p  >  a,  the  Fourier  series  expansion  for  E 

inC 

consists  of  the  Fourier  series  expansion  of  E and  a  sum  of  outward- 

xnc 

traveling  waves.  Here,  E  is  the  z  component  of  the  incident  electric 

z 

xnc 

field.  The  x  coordinate  that  appears  in  expression  (120)  for  jE  is 
expressed  in  terms  of  o  and  <J)  as 


x  *  a  cos  i>  +0  cos  P 
o 

Substituting  (B-3)  into  (120),  and  then  taking  the  z  component  of  the 

resulting  equation,  we  obtain 

,  -jka  cos  d>  .. 

_inc  J  o  -jko  cos  p 

E  =  e  e 

z 

If  [9,  eq.  (5-101)]  is  used  to  expand  the  second  exponential  in  (B-4), 
then  (B-4)  becomes 


(B-j 


(B-* 


Einc  -  e 
z 


-jka  cos  2  30 


O  v  ,-n  T  n  ,\  ln® 
l  j  Jn(ko)eJ 


(B-5 


n=-=° 

where  J  is  the  Bessel  function  of  the  first  kind  of  order  n. 


Since  E  is  finite  at  p  »  0,  the  Fourier  series  expansion  for  E  in 
z  z 

the  region  for  which  p  <  a  can  have  only  Bessel  functions  of  the  first 

kind.  Moreover,  the  Fourier  series  expansion  for  Ez  in  the  region  for 

Xnc 

which  o  >  a  consists  of  Ez  of  (B-5)  and  a  sum  of  outward  traveling  waves. 
Therefore, 


(“  -jka  cos  <p  00 


0  l  A  J  (ko)ejnd> 
n  n 


p  <  a 


(B-6 


-jka  cos  <$>  ®  .. 

°  E  (j  Jn(kp)  +  BnH^2)(kp))ejn4>, 

n«-« 


e 


o  >  a 


where  H  '  is  the  Hankel  function  of  the  second  kind  of  order  n,  and  A 
n  n 

and  B  are  unknown  complex  constants.  Requiring  E  of  (B-6)  to  be  con- 
n  2 

tinuous  at  0=a,  we  obtain 


B  = 
n 


(A  -  j“n)J  (ka) 
n _ n 

H(2)(ka) 

n 


(B-7) 


Substitution  of  (B-7)  into  (B-6)  gives 


-jka  cos  j>  00  j.. 

.  0  l  A  J  (kp)eJ"* 

L  n  n 
n=-°° 


E  =< 
2 


I  -jka  cos  4> 


l  (j-nJn(ko)  +-a 


,  0  £  a 

(B-8) 

(A„-rn>  J  0«>  (2)  jn4) 

— -  HU;(kp))ejnQ,  p  >  a 

(2)  v  n  — 


Hv  '(ka) 
n 


Requiring  the  upper  right-hand  side  of  (B-8)  at  p=a  to  be  equal  to 
the  right-hand  side  of  (B-l),  we  obtain 


-jka  cos  <j> 
e 


CO 


I 


A  J  (ka)e^m<^  = 
m  m 


(B-9) 

other  values  of  6 


The  new  summation  index  m  in  (B~9)  was  chosen  for  convenience.  Multiplying 
both  sides  of  (B-9)  by  e-'^  and  integrating  from  0  to  2tt  with  respect  to 
<p,  we  obtain 


-jka  cos  j) 

2tt  A  J  (ka)e 
n  n 


V 


180°+<i> 

o 

l80°-4> 

o 


Ji  -  d* 

^o 

(B-10) 


If  the  quantity  in  parentheses  on  the  right-hand  side  of  (B-10)  is  called 


x,  then  (B-10)  reduces  to 
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•  -  *•  _  »  L  «•  i."W  .  ".  ■*_  “-  »v~ -  -^  *.’.  <  ”\r  . 


^ar  ^yi 

fc&&2 


According  to  [12,  Formula  3.752  (2.)], 
1 


f- 


x  cos  (n<{>  x)dx  =  tt 


V^V 


i  n  I  4>, 


(B-12) 


-1 


Equation  (B-12)  is  written  with  the  understanding  that  J,  (|n|<J>  )/(!n|<b  )  is 

1  o  o 

to  be  replaced  by  1/2  when  n  is  zero.  Substitution  of  (B-12)  into  (B-ll) 
leads  to 


A  J  (ka)e 
n  n 


-jka  cos 


V<J>  (-1)  J.  ( J  n )  4>  ) 


(- 


1 x  1  ■  ^o 


n|<j> 


) 


(B-13) 


Substituting  (B-13)  into  (B-8),  we  obtain 


_ !!»  ?  MfL  -Ji(|n|V-  , 

?  ^  J  (ka)  n | $  n 


n=-°o  n 


(ko)e^n<^  ,  o  <  a 


E  = 
z 


J 

\ 


V<t> 

— 2 


°  I 


.  ,  .n  J_  (  j n  j 4)  )  ,ON 

&rt; ( ■  ini^-0")  “iS  *  + 


(B-14) 


n=-°°  H  (ka) 
n 


L 


-jka  cos  $  00  J  (ka) 

or  •  -n .  T  ,  n 
i  I  j  (Jn(kp)  -  — 

n=-oo 


H(2) (ka) 
n 


H^2)(kD))ejn(b,p  >  a 


With  a  view  toward  enforcing  (B-2) ,  we  call  the  d)  component  of  the 


magnetic  field  and  express  it  as  [9,  eq.  (5-18)] 


.  3E 

H  =  =1  — 5- 

f  kn  9p 


(B-15) 


The  limits  of  (B-15)  as  p  approaches  a  from  below  and  above  are 


H"  =  -  J-  lim 
*  kn  p^a”  ap 


(B-16) 


and 


,  .  SE 

Ht  *  -  J-  iim  — 5. 

*  kn  p-va+  30 


(B-17) 


hMk 

.-Vv' 

•:-:nva 


av- 


‘W 


i  *  A 
*  .  ■  ,  * 


r<u- 


-.VV 


v.va). 

c  --  •• 


,v-y* 


W 


sn»> 


.s  , 


.  *  »  •  »  .  .  .  • 


where  lim  and  lim  denote  the  limits  as  P  approaches  a  from  below  and 
fH-a"  0->-a+ 

respectively.  Substitution  of  (B-14)  into  (B-16)  and  (B-17)  gives 


and 


_  jV4> 

"  ~~TrT  l 


(-l)n  ^Ji  C '  n  1 0p) 


4>  2n  4  j  (ka) 

n=^=°  n 


M<J> 


■)  J'(ka)e 
n 


Jn4> 


n  _ 

H(^  2n  n--«  H^2)(ka)  j n | 4>q 


-jka  cos  <j>  00  J  (ka)  , 


I  j“n(J’(ka)  -  -Ssr -  (ka))eJ 

n**-00  n  '  (ka) 

n 


The  identity  [9,  eqs.  (D-12)  and  (D-17)] 


V*)  “n^00  ‘  Jn(*)  an2)<lt)  ’  ”  TO 


reduces  the  difference  of  (B-18)  and  (B-19)  to 


,  Vcf>  <”  J.  ( !  n  1 4>  ) 

<  - »;  -  ^  i  w>v — *->  <• 

n=-°° 


J  (ka)  H(2)(ka) 
n  n 


)  P"*  + 


„  -jka  cos  $  00  ,-n  jn<j> 

2  °  S  J  e 

^kae  naL,  H(2)(ka) 

n 


Substitution  of  (B-21)  into  (B-2)  and  use  of  the  identity 
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